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Abstract: The globular clusters are probably good targets for dark matter (DM) searches 
in 7-rays due to the possible adiabatic contraction of DM by baryons. In this work we 
analyse the three-year data collected by Fermi Large Area Telescope of globular clusters 
NGC 6388 and M 15 to search for possible DM signals. For NGC 6388 the detection of 
7-ray emission was reported by Fermi collaboration, which is consistent with the emission 
of a population of millisecond pulsars. The spectral shape of NGC 6388 is also shown 
to be consistent with a DM contribution if assuming the annihilation final state is hh. 
No significant 7-ray emission from M 15 is observed. We give the upper limits of DM 
contribution to 7-ray emission in both NGC 6388 and M 15, for annihilation final states 
66, W^"*"!^", T^T~ and monochromatic line. The constraints are stronger than that 

derived from observation of dwarf galaxies by Fermi. 
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1. Introduction 

A standard model of cosmology is developed, in which the universe consists of 4% ordinary 
baryonic matter, ~ 23% dark matter (DM), ~ 73% dark energy, and a tiny abundance 
of relic neutrinos Q. The nature of DM particle remains a mystery. One of the leading 
candidates is the weakly interacting massive particle (WIMP), which is predicted in several 
models, such as neutralino in supersymmetry model (see the reviews ||2|, ^, ^). In this kind 
of models, the mass and interaction strength of DM particles can produce the correct relic 
density of DM if the WIMPs are thermally "freeze-out" , which is called "WIMP miracle" . 

If DM particles annihilate or decay into standard model particles, they can be detected 
indirectly from the cosmic ray (CR) radiation. Among many kinds of CR particles, 7-rays 
are the best probe due to their simple propagation. Fermi gamma-ray telescope, which 
was launched in 2008, has surveyed the 7-ray sky with very high resolution and sensitivity 
for more than three years. Nearly 2000 sources as well as the diffuse 7-ray emission were 
detected by Fermi Large Area Telescope {Fermi-L AT) ^, 0, ^] . The analysis of the Fermi- 
LAT data in the Galactic center region did see some excesses with respect to the background 
model IC], however, there is no strong indication of signals from DM annihilation or 
decay^. The constraints on DM model parameters can be derived according to the non- 



detection of DM signals from e.g., dwarf galaxies I|l2|, |13|, |15|, |16|, galaxy clusters ||17 

H, n, 1^, n 



and the diffuse 7-rays |2J, 



25, 29, pq, 31 



Si- 



Due to the very weak interactions of DM particles, it is important to investigate the 
sites with high DM density when searching for DM annihilation signals. The proposed 



^See also the argument of possible DM explanation of the 7-ray haze/bubble 
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good candidates include the Galactic center, dwarf galaxies, Galactic subhalos and cluster 
of galaxies. The Milky Way globular clusters (GCs), defined as spherical ensemble of stars 
that orbits the Galaxy as satellites, are also potentially good targets for indirect detection 
of DM. The formation of GCs remains a poorly understood problem. There are generally 
two scenarios to describe the formation of GCs. The primordial formation scenario suggests 
that GCs were formed in cosmological DM minihalos before the formation of galaxies p^ . 
The other way to form GCs might be the star-forming events such as the merger of galaxies. 
There was evidence to show that metal-poor GCs might have a cosmological origin and 
metal-rich GCs might form in the galaxies |35]. If the GCs were formed in cosmological 
DM minihalos, they would experience the adiabatic contraction (AC) due to the infall of 
baryons during the evolution of GCs and leave a high density spike of DM. GCs are not 
usually discussed for DM detection due to the poor knowledge about their origin and the 
observational fact that there is in general no significant amount of DM in vicinity of GCs 
1 36, 38 1 . However, there is possibility that the high density spike of DM due to the AC 
process may still play an important role for the annihilation signals. The previous works 
to search for or constrain DM models with 7-rays from GCs include |3^, |40| . 

Recently the atmospheric Cherenkov telescope array High Energy Stereoscopic System 
(H.E.S.S.) had investigated two GCs NGC 6388 and M 15 to search for possible DM signals 
||4l[] . No 7-ray signal was detected by H.E.S.S. and strong constraints on the DM model 
parameters were given. In this work, we use the three- year data of Fermi-hAT to study the 
7-ray emission from DM annihilation in these two GCs. Detections of 7-ray emission from 
some GCs with Fermi-L AT were reported [42, 43, 44, including NGC 6388 studied here. 
For M 15 there is no detection yet. In this work we will focus on the possible DM component 
of the 7-ray emission, if any, from the GCs. The upper limits of DM contribution will be 
derived and the constraints on DM model parameters will be presented. 



2. Gamma-rays from DM annihilation in globular clusters 



M 15 is a metal-poor GC which favors a cosmological origin of it |46|| . For NGC 6388, 
there is strong evidence to show the existence of an intermediate mass black hole (IMBH) 



with mass ~ 6 x 10^ Mq |47], which also suggests a cosmological origin even though the 
metallicity is relatively high [|^. Therefore we have good motivation to search for the 
possible DM annihilation signal from these two GCs. The estimated stellar masses of NGC 
6388 and M 15 are 10^ and 5 x 10^ Mq, with distances 11.5 and 10.0 kpc respectively 
1 47, 40 1 . Other parameters of them can be found in Table 2 of Ref. |41|. 



2.1 DM density distribution 

For the purpose of this work, these two GCs are assumed to form in the cosmological 
context, which were DM dominated in the primordial stage, before reionization and the 
galaxy formation . The AC process of baryons to form the GC is expected to pull DM 



into the center and results in a high density core of DM |48]. After the AC process the 



heating effect of DM due to scattering with baryons will tend to sweep out the high density 
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DM core, leaving a constant density ||49[]. The IMBH, if exists, may further modify the 
density profile through adiabatic accretion [ pC[| . 

The modelling of the GC DM halo can be divided into three steps. The first step is 
the AC process of the dark halo during the collapse of the core of GC. Supposing that the 
DM particles travel on circular orbits, the enclose mass distribution of DM M{r) can be 
calculated with the follow equation 

[MDM,i(ri) + Mb^i{ri)]ri = [M^Mjirf) + iV4j(rj)]ry, (2.1) 

where the subscript i{f) denotes the initial (final) mass distribution of baryon or DM. The 
initial mass of the minihalo is assumed to be 10'' Mq, with Navarro- Frenk- White (NFW, 
|51]) density profile for both the DM and baryon distributions^. The mass fraction of 



baryons is adopted to be 20%. For the convenience of comparison, these adoptions are the 



same as that in Ref. |41|. We should keep in mind that these parameters may have large 
uncertainties and the quantitative results of this work may also suffer from uncertainties. 
Given the final baryon distribution, which can be derived according to the observed 



surface density distribution of the GC^, one can get the DM density profile after AC |48|. 
The final baryon density for NGC 6388 is taken from [^], which was computed using the 
surface density profile given in ||4^. For M 15 the final baryon density is taken from |[55[| . 

The second step is to take into account the smoothing effect due to baryon heating 
after AC process. For the convenience of discussion, we employ the relaxation time Tr 



defined as |56] 



3.4 X 10^ / Urms \^ f m- \ ^ f n ^ 



where frms is the velocity dispersion of stars, m is the typical stellar mass in the GC, n 
is the stellar number density, and In A is the Coulomb logarithm. Tr is estimated to be 
~ 7 X 10^ yr in the central region of M 15, and ~ 8 x 10^ yr for central NGC 6388 |l|. 
The relaxation time is an increasing function of the distance to the center. The DM will 
be heated up due to the scattering with the stars, which will lead to the dissipation of the 
DM core [|^. The scattering time scale is comparable to T^. The heating radius, rheat, is 
then defined with T,. (rheat) = ^age, where tage is the age of the Universe. Therefore at small 
radius r < rheat ) the relaxation time is shorter than the age of the Universe and the heating 
effect on DM is important. At large radius r > rheat the DM distribution is unaffected by 
heating. The heating radii are estimated to be about 5 pc and 4 pc for M 15 and NGC 
6388 respectively. Roughly speaking we have the DM density distribution with baryonic 
heating 

. . / ^0, r< rheat, , . 

[ ^ PAc(r-hcat)' ^ > ^I'^'^t' 

where pAci''') is the DM density profile after AC, which can be solved with Eq. (2.1), and 
po is the density at r = rheat- 



^Note that for such minihalos the density profile might be smoother however, as shown in jHi 

the initial density profile does not affect significantly the final DM profile after AC process. 
^See, e.g., http://www.physics.mcmaster.ca/ harris/mwgc.dat 
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The third step is to consider the effect of the IMBH, if exists. The AC profile of DM 
will not be significantly affected by the IMBH due to its small mass compared with the 
total baryon mass. However, the following adiabatic accetion of IMBH will modify the DM 
density profile after dynamic heating. The radius within which the IMBH is gravitational 
dominant is defined by M{< r^) = /j(r)d^r = 2Mimbh- Then the DM distribution can 
be expressed in three regions 



Po{r/rh) r < rh, 

Po, rh<r < rheat, (2.4) 

k PAClf-hoat) ' ^^^^ 



The inner most density profile (oc r ^/^) corresponds to the collisionally regenerated struc- 
tures ("crest") of DM due to the joint evolution of baryons and DM in the enviroment of 
a central black hole ||]. For NGC 6388, Mimbh ~ 6 x 10^ Mq, r,, ~ 0.4 pc is found with 
the observational baryon density. For M 15 there is no consensus of the existence of IMBH 
1 59, pO], and we employ Eq. (2.3) to describe the DM density profile of M 15 without 



considering the possible IMBH. 
2.2 Astrophysical J-factor 

For Majorana fermion DM particles, the 7-ray flux from DM annihilation can be written 
as 

$(AJ7, E,) = j-x^^x J{An)A^, (2.5) 

where niy. and {av) are the mass and velocity weighted thermal average annihilation cross 
section of DM particles, dNy/dE^^ is the 7-ray spectrum for one annihilation. The as- 
trophysical factor (J) is the integral of the density square along the line of sight (LOS) 
averaged over the solid angle Aft 

J(Af]) = -!-/" dn [ dlp\r{l)). (2.6) 
For H.E.S.S. observations, the integral solid angle is Ail = 5 x 10^^ sr, which corre- 



sponds to a cone with half angle 0.07° [41|. Since the resolution angle of Fermi-LAT in GeV 



range is much larger (> 0.5°, |61[]), we need to enlarge the integral solid angle. The tidal 



radii of both GCs are about 30 pc, and the distances are about 10 kpc [41j. The opening 
angles of these two GCs are ~ 0.17°. Therefore they can be regarded as point sources for 
Fermz-LAT and we integrate all the DM contribution to the tidal radius to calculate the 
J-factor. It is found that the final J x AJl is about 7.8 x 10^^ (3.4 x lO^O) GeV^ cm'^ for 
M 15 (NGC 6388), which is larger by ^ 10% (0.1%) compared with that within 0.07° cone 
as adopted by H.E.S.S.. The J-factor of NGC 6388 is larger than M 15 is mainly due to 
the difference of the density profiles after AC which depends on the final baryon density 
profiles of the GCs, and the heating effect. Compared with the dwarf galaxies as given in 
p5[], the J- factors of GCs are generally larger, which is also due to the AC process. 
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3. Fermi-'LAT data analysis 



The Fermi-LAT data^ used in this analysis are the new "Pass 7" data recorded between 
4 August 2008 and 2 September 2011. Photons with Event Class "Source" (evclass=2) 
and zenith angle within 100° are selected. The energy range of events is cut from 200 
MeV to 300 GeV, and the radius of region-of- interest (ROI) is adopted to be 6°. We 
use the LAT Scientific Tools v9r23pl to do this analysis. The unbinned likelihood anal- 
ysis method is adopted. The instrument response function used is "P7S0URCE_V6" . For 
the diffuse background, we use the Galactic diffuse model gal_2yearp7v6_v0 . f its and 
the isotropic background spectrum iso_p7v6source.txt provided by the Fermi Science 
Support Center^. 

3.1 NGC 6338 



The detection of 7-ray emission from NGC 6388 was reported in |43|, with Test Statistic |62| 
value TS= 86.6. The spectral energy distribution can be fitted using power-law function 

with an exponential cutoff ex.p{—E /Ec), which is expected for the emission from a 

1-0.7 
-0.5 



population of millisecond pulsars (MSP). The best-fitting parameters are T = l.l^^'l^ and 



Ec = l-8lo:? GeV H]. Using the hkelihood tool gtlike we re-do the spectral analysis 
with more data. The source model XML file is generated using the user contributed tool 
make2FGLxml . py^ based on the 2FGL source catalog ||5|. The spectrum of NGC 6388 is 
also modelled with E~^ exp(— E'/i^c)- By setting all the source parameters within the ROI 
free, the best-fitting parameters for NGC 6388 are T = 1.21 ± 0.17 and E^ = 1.82 ± 0.35 



GeV, with a TS value 596. The fitting parameters are consistent with that given in [43|. 

To derive the spectral energy distribution (SED) of NGC 6388, we divide the data 
into different energy bins, and use gtlike tool to fit the parameters for each bin. Two 
methods are adopted in the fit. We first fix the parameters of all other sources and the 
normalizations of diffuse backgrounds derived above in the global fit, leaving only the 
normalization parameters of NGC 6388 and the very bright pulsar PSR J1709-4429 free. 
The spectral parameters of NGC 6388 and PSR J1709-4429 are also fixed to be the best 
fitting values. Because the energy bin is relatively narrow the precise values of the spectral 
parameters have little effect on the final results |^3[. The results are shown by the filled 
circles in Figure [l|. The solid line shows the fitting curve with spectrum E~^ ex.p{—E/Ec), 
as a represent of MSP-type emission. It shows good agreement between the global fit and 
the individual fit for different energy bins. Then we relax the normalization parameters 
of all the sources and the diffuse backgrounds, and re-do the fit. The results are shown 
by the empty circles in Figure |. We can see that there are some differences between the 
results of these two methods, which might be originated from the complexity of the diffuse 
background models. 

We then add a DM component at the position of NGC 6388 to search for possible DM 
contribution to the emission. Different annihilation final states are investigated, including 



* http : //fermi. gsfc. nasa. gov / ssc /data 

^http:/ /fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels. html 
®http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/user/ 
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Figure 1: SED of NGC 6388. 



bb, W~^W , fi'^fi , T^T and 7-ray line. We use PYTHIA simulation tool to calculate the 
7-ray yield spectrum for each final state [33|. A series values of DM mass from 10 GeV to 



10 TeV are considered. For monochromatic 7-ray line we use Gaussian function to model 
its spectral shape. The energies of photons from 300 MeV to 200 GeV are searched. The 
width of Gaussian function as a function of photon energy is adopted to be the energy 
resolution of Fermi-L AT, which is about 8% — 13% (AE/E) in this energy range |l6l| . 
Given the spectrum shape of DM contribution, we use the python likelihood tooF to fit 
the normalization parameter and derive the fiux upper limits. 

It is found that the detected 7-ray spectrum of NGC 6388 can also be fitted with a DM 
component, with mass ~ 25 GeV and annihilation final state bb, as shown by the dashed 
line in Figure |l|. Note that the recent analysis of the Fermi-LAT data in the Galactic 
center region also showed possible additional emission compatible with DM contribution 
with ~ 30 GeV for bb annihilation final state [64|. 



It is well motivated that the 7-ray emission from GCs may come from the MSPs, 
therefore we do not claim the DM origin of the 7-rays of NGC 6388. In any case we can 
instead set an upper limit of the contribution from DM annihilation. Here the upper limits 
are derived for different final states and diff'erent mass of DM particle individually. We 
use two ways to fit the data. The first one (Fit 1) is to fix all the parameters of sources 
derived in the above global fit. The free parameters are the normalizations of the diffuse 
backgrounds and the DM component. The other way (Fit 2) is to leave the normalizations 
of sources in the ROI and the diffuse backgrounds as well as the normalization of the DM 
component free. The spectral parameters of the sources are fixed to be the best fitting 
values in the global fit. Fit 1 corresponds to a more stringent constraint, based on the 
assumption that the observed 7-rays come from astrophysical sources. Fit 2 gives a weaker 
but more conservative constraint. The 95% confidence level (C.L.) upper limits of the 
> 100 MeV fluxes from DM annihilation for NGC 6388 are shown in the left panel of 



''http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/python_tutorial. html 
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Figure The thick and thin hnes correspond to Fit 1 and 2 respectively. The same way 
is apphed to the Hne analysis. The upper limits of line emission are shown in Figure ^. 




Figure 2: Derived 95% C.L. upper limits of the WIMP annihilation contribution to the > 100 
MeV 7-ray fluxes as functions of WIMP mass m^, for GCs NGC 6388 (left) and M 15 (right). 




Figure 3: 95% C.L. upper limits of the monochromatic 7-ray emission as functions of 7-ray energy. 



3.2 M 15 

There was no firm detection of emission from M 15 in the previous analysis. In [^ ] M 15 
was reported to have a very weak signal with TS= 5.4. In our analysis with more data, we 
find a relatively higher TS value ~ 12, for both power-law and power-law -|- exponential 
cutoff models. We can derive the upper limits on DM annihilation to 7-rays, similar as the 
analysis of NGC 6388. Because there is no detection of 7-rays, only the method Fit 2 is 
adopted for M 15. There are three other sources in the ROI of M 15, 2FGL J2115. 4-^1213, 
2FGL J2112.5-F0818 and 2FGL J2147.3-h0930. The free parameters include the spectral 
parameters of M 15 (power-law model for a possible MSP contribution) and these three 
sources, the normalizations of the diffuse backgrounds and the DM component of M 15. 
The results are shown in Figures I and |. 



-7- 



4. Constraints on DM models 



Integrating Eq. ( |2.5| ) above 100 MeV we can easily get the upper limits of the DM anni- 
hilation cross section using the flux upper limits. The results are presented in Figures ^ - 
|6| respectively. 

Figure^ shows the constraints on m^ — {av) for hh (left) and W~^W~ (right) final states. 
For comparison we also show the results from the combined analysis of Fermi observations 
of 10 dwarf galaxies [|l5| and that given by H.E.S.S. observation of these two GCs 
||4l|] ). It is shown that the constrains given in this work are generally stronger than that 
for dwarf galaxies, at least for DM mass > 50 GeV. Because of the high density spike 
from the AC process, the distribution of DM is more concentrated than the initial NFW 
profile. Although the heating effect from stars will smooth out the central density spike, 
the J- factors of GCs are still higher than that of dwarf galaxies (e.g., estimated with NFW 
profiles), and the constraints on DM models are stronger accordingly. For bb final state, the 
DM annihilation induced 7-ray spectrum is similar with the observed data of NGC 6388, 
therefore the constraint is a bit weaker when < 50 GeV. However, for the method Fit 
1 the constraint is always stronger than that for dwarf galaxies. H.E.S.S. constraints are 
more effective for massive DM {my. > 2 TeV). 

Also shown in Figure ^ are the theoretically expected neutralino annihilation cross 
sections (multiplied with the branching ratios to bb and WW + ZZ) in the Minimum 
Supersymmetric Standard Model (MSSM). In the right panel we sum the model predicted 
cross sections to W~^W~ and ZZ chanels together due to the similarity of 7-ray spectra 
from these two channels. We utilize numerical code micrOMEGAs []65| , |66| to perform 
a random scan in the 7-dimensional parameter space at the electroweak scale. These 
parameters include the CP-odd Higgs mass m^, the Higgs mixing mass parameter the 
wino mass parameter M2, the sfermion mass parameter mj, the ratio of two Higgs vacuum 
expectation values tan /?, the trilinear parameters of the third family squark and Ai. The 
other trilinear parameters are set to zero. We also impose the assumptions that the gaugino 
mass parameters are related by Mi : M2 : M3 = ai : a2 '■ o-d, for grand unification, where 
the ai are the coupling constants of three standard model gauge groups, and Mi : M2 = 
|tan2 6',^. The ranges of the parameters are taken as follows: 50 GeV< |/i|,M2 <10 TeV, 
100 GeV < mA,mj: <1 TeV, 1< tan(3 <60, -bnij: < At,Ab < Bmj and sign(//) = ±1. 
Several constraints from accelerator experiments and DM detection are implemented in 
our numerical scan. We set the limit for p parameter as p — 1 < 2.2 x 10~^|^^. Some 
important flavor physics constraints include: Br{B — > ^^7) = (3.55 it 0.24) x 10~^ pS)] , 
Br{Bs /x+/i-) = (0 ± 1.4) X 10-® ||, 0], Br{Bu Tv)/Br{Bu -f tv)sm = 1-28 ± 0.38 
||6^ . Here we only require the supersymmetric contributions to satisfy these constraints 
at 3a" level, and adopt a very conservative bound for muon anomalous magnetic moment 
|7l|] as -11.4 X 10"^° < 5a^ < 9.4 x 10"^ |6^. The mass bound of standard model like 
Higgs rufi > 114 GeV, limits on the masses of light charge sparticle from LEP (for details, 
see [S5, |6^), and DM direct detection constrains from XENONIOO [72| are also taken into 
account. 

In Figure Q the squares are for DM models which can give the right relic density Q 
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if DM is thermally produced in the early Universe (labelled as "MSSM-thermal"). The 
triangles are the cases with thermal relic density not higher than the measured value, the 
correct DM abundance in these models could be produced via some non-thermal mech- 
anisms. We can see that for bb final state the current constraint can reach the natural 
scale for thermally produced DM with mass of 0(10) GeV. In the MSSM model, large 
neutralino annihilation cross section to bb final states could arise from the resonance effect 
with niA ~ 2mDM- Some non-thermal models with this feature have been excluded by 
our constraints. For W~^W~ channel the constraint is a bit weaker but also close to the 
natural scale. In the MSSM scenario, large higgsino or wino component (depending on the 
relations between the three parameters Mi, M2 and fi) in the neutralino would enhance 
DM annihilation cross section to gauge bosons significantly. If the neutralino mass lies in 
the range (80, 300) GeV, many non-thermal models would be also stringently constrained 
by our results. 
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Figure 4: Constraints on DM mass vs. annihilation cross section, for bb (left) and W^W^ (right) 
final states. Points are a random scan of the MSSM parameter space taking into account the current 
constraints from accelerator data. Magenta dotted lines in the right panel are the constraints got 
by H.E.S.S. observations of NGC 6388 (lower) and M 15 (upper). 




Figure 5: Same as Fig. ^ but for leptonic final states (left) and r+r (right) respectively. 



Motivated by the recent observations of the CR positron/electron excesses at PAMELA, 
ATIC and Fermi-LAT [||, 0, ^, ^ and the non-excess of antiprotons p|, the 
leptonic DM models are proposed to explain the data (e.g. 



8C, 31, 82, B3, 84 



We 
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Figure 6: Same as Fig. ||but for DM annihilatfon into 77 (left) and 7^° (right) final states. 



also study the constraints on the leptonic final states jU and t~^t~ by DM annihilation 
■which might be responsible for the excesses. Here the inverse Compton scattering 7-rays 
generated by the decaying products from muons or tauons are not taken into account. 
As illustrated in ||l^ , considering the inverse Compton 7-rays the constraints on DM model 
parameters will be stronger, however, depending on the uncertainties of the diffusion pro- 
cess of e^. Thus the results given here should be conservative. Constraints on — (av) 
parameter space are shown in Figure ^. The contours show the favored parameter region 
to fit the CR data [33|. It is shown that the models to explain the excesses should 
be excluded by the Fermi-LAT data about GCs. 

Finally we study the constraints on possible monochromatic 7-ray line emission from 
6-g-, XX ~^ 77 XX ~^ 7^ of DM annihilation. No significant line emission is found 
in the data. The upper limits of line emission are derived. The constraints on cross 
sections to 77 and 7Z are given in Figure |^. Note we have E-y = niy. for XX ~^ 77) ^-iid 



m^(l 



m^/Am'r) for xX ~^ 1^- The results derived with Fermz-LAT data including 



the Galactic center region by Fermi collaboration (NFW profile, [£6|) and Vertonger & 
Weniger [37| are also shown for comparison. Our constraints are a bit weaker than the 
results in these two works. We think it is reasonable because their analysis regions are 
much larger and include the Galactic center region, which will give a higher J-factor of 
DM annihilation. 



5. Conclusions and Discussions 



The GCs are thought to form in the cosmological context with AC process at the beginning 
which pulls DM into the halo center and results in a very high annihilation luminosity. Thus 
search for 7-rays from GCs may be effective to probe the particle nature of DM. In this 
work we analyze the Fermi-LAT three-year data (Pass 7) of GCs NGC 6638 and M 15 
and constrain the DM annihilation models. A clear detection of 7-ray emission from NGC 
6388 is found, with TS value ^ 600. The spectrum of NGC 6388 can be well fitted with a 
power-law + exponential cutoff function, which is expected for the emission of a population 
of MSPs. We find that a DM scenario with ttIv ^ 25 GeV and hh final state can also fit the 
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SED. For M 15 no significant 7-ray emission if found (the spectral fit indicates a potential 
source with TS ^ 12). 

Assuming there is an additional spectral component from DM annihilation of these two 
GCs, we derive the upper limits of the DM component for different DM masses (lOGeV — 
lOTeV) and annihilation final states {bb, W~^W^ , j-i^fJ-^, t^t^, 77, jZ) (Figures § and |3|) . 
The constraints on the DM annihilation cross section are given (Figures ^[^. Except for 
the line emissions, the constraints are stronger than that derived according to the Fermi- 
LAT observations of dwarf galaxies. For DM mass smaller than TeV our constraints are 
also stronger than that given by H.E.S.S. observations of the same GCs. For bb and W^W~ 
final states which are generally expected from supersymmetric DM model, the constraints 
can reach the natural scale with which DM is thermally produced. Especially the leptonic 
annihilation models to explain the CR excesses can be excluded by the current analysis. 

However, the uncertainties of the present analysis, for example the properties of the 
hypothetical DM halo and the origin and evolution of the GCs, are far from clear. The 
GCs were assumed to be formed in the cosmological context, and the DM density profiles 
in the GCs are modelled taking into account the most probable astrophysical processes, 
e.g., the AC by baryons, adiabatic growth of an IMBH and the scattering by stars. Future 
studies on the observations and modelings of the DM distribution in the GCs are necessary 
to improve the current work. 
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